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The statistics of Lagrangian pair dispersion in a homogeneous isotropic flow is investigated 
by means of direct numerical simulations. The focus is on deviations from Richardson eddy- 
diffusivity model and in particular on the strong fluctuations experienced by tracers. Evidence 
is obtained that the distribution of distances attains an almost self-similar regime character- 
ized by a very weak intermittency. The timescale of convergence to this behavior is found to 
be given by the kinetic energy dissipation time measured at the scale of the initial scpara- 
I tion. Conversely the velocity differences between tracers are displaying a strongly anomalous 

behavior whose scaling properties are very close to that of Lagrangian structure functions. 
' O ' These violent fluctuations are interpreted geometrically and are shown to be responsible for 

a long-term memory of the initial separation. Despite this strong intermittency, it is found 
that the mixed moment defined by the ratio between the cube of the longitudinal velocity 
difference and the distance attains a statistically stationary regime on very short timescales. 
These results are brought together to address the question of violent events in the distribution 
of distances. It is found that distances much larger than the average are reached by pairs that 
have always separated faster since the initial time. They contribute a stretched exponential 
■ behavior in the tail of the inter-tracer distance probability distribution. The tail approaches 

a pure exponential at large times, contradicting Richardson diffusive approach. At the same 
' time, the distance distribution displays a time-dependent power-law behavior at very small 

values, which is interpreted in terms of fractal geometry. It is argued and demonstrated nu- 
merically that the exponent converges to one at large time, again in conffict with Richardson's 
distribution. 
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1. Introduction 

It is known since Taylor's seminal work that tracers transported by a turbulent 
flow approach a diffuse behavior on time scales much longer than the Lagrangian 
correlation time of the flow. These ideas are now commonly used in applications, 
as for instance in air quality control, to model effective mixing properties in terms 
of an eddy diffusivity. Such models give a good handle on long-term averages and 
are successfully used to determine, for instance, possible health hazards linked to a 
long exposure downstream a pollutant source. However they are unable to capture 
strong local fluctuations stemming from the complex structure of the turbulent 
flow. Such events cannot be directly predicted from the averaged concentration 
field as they relate to higher-order moments. Accessing these fluctuations is crucial 
in order to quantify for instance the likeliness of finding a local concentration 
exceeding a high threshold. 

Second-order statistics, such as the variance of a transported concentration field 
and more generally the spatial correlation of a passive scalar, are statistically re- 
lated to the relative motion of tracers (see, e.g., The problem consists then in 
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investigating the time evolution of the separation R(t) = Xi{t) — X2{t) between 
two Lagrangian trajectories. In turbulence, the distance |-R| follows Richardson- 
Obukhov super diffusive law ~ et^, where e is the mean rate of kinetic 

energy dissipation. The long-term behavior is thus becoming independent of the 
initial separation |-R(0)| = ro, whence the designation of explosive pair separa- 
tion. In his original arguments to derive this law, Richardson assumed that 
pair separation is a diffusion process with a scale dependent diffusivity K. Then 
the transition probability density p(r,t\rQ,0) that two tracers are at a distance 
= r knowing that they were initially separated by |-R(0)| = ro satisfies the 
Fokker-Planck equation 



dtp = -^dr [r'^K{r) drp] 



(1) 



The atmospheric measurements of Richardson, later refined by Obukhov [4] in view 
of Kolmogorov 1941 phenomenology, led to assume that for separations r within 
the inertial range of turbulence K{r) oc e^/'^r^/^. This implies that, at large times, 
oc et^ and the transition probability reads 



p{r,t\ro,0) oc 



|i?(t)|2)3/2 



exp 



^^2/3 



|i?(t)|2)l/3j 



(2) 



where ^ is a positive constant. The considerations leading to the Fokker-Planck 
equation ([1]) rely on the use of a central-limit theorem for R{t) and thus on the 
hypothesis that the velocity difference between the two tracers is correlated over 
timescales much smaller than those of interest. As noticed for instance in such 
an hypothesis can hardly be invoked. It is indeed known in turbulence that eddies 
of size r are correlated over a time of the order of their turnover time ~ ^-1/3^,2/3^ 
Hence for separations that grow like r ~ eV2i3/2 

one has r,. ~ i , so that the dom- 
inant flow structures in the separation dynamics are in principle correlated over 
timescales of the order of the observation time. Despite such shortcomings, the 
diffusive approach proposed by Richardson and in particular the explicit form ([2]) 
for the transition probability density have proven being relevant in some asymp- 
totics Also, much work on relative dispersion has used it as a basis. For 

instance, improvements of ([2]) were proposed using modified versions of the eddy 
diffusivity K{r), adding a time dependence [ill . Il2i] or, more recently, including 
finite Reynolds number effects [l3]. All of these improvements strongly alter the 
functional form of the large-r tail of the transition probability density. Neverthe- 
less, the physical mechanisms leading to the de-correlation of velocity differences 
and to models based on eddy diffusivity are still unclear and many questions re- 
main open. For instance, the speed of convergence to the Richardson-Obukhov law 
and the form of subleading terms are still not known. 

The first work dealing with the way pair separations converge to a superdiffusive 
behavior is due to Batchelor [ij]. He argued that the explosive law is preceded by 
a ballistic phase during which the tracers keep their initial velocity and separate lin- 
early in time, i.e. R{t) — ro ~ tSrgU ~ i(ero)"'^/^, up to a time Trg ~ e'^^^r^^^, which 
is equal to the eddy turnover time associated to the initial separation. This first 
regime where velocity remains strongly correlated can clearly not be described by 
eddy-diffusivity approaches. Various stochastic models have been designed to catch 
the two regimes. Most of them are based on the observation that the acceleration 
difference between the two tracers is shortly correlated and do not assume that ve- 
locity differences get uncorrelated. The pair separation and the velocity difference 
can then be approximated as coupled Markovian diffusive processes (see [15|, ll6|] 
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for reviews). The usual path for designing such models requires imposing some 
constraints on the drift and the diffusion terms. Thomson [17l | argued that they 
should satisfy the well-mixed condition: when averaging over uniformly separated 
pairs inside the whole inertial range, the statistics of velocity differences between 
the two tracers has to recover Eulerian two-point statistics. Definitively devising 
an admissible model requires an input from Eulerian statistics [18]. It is known in 
turbulence that the distribution of velocity differences is neither self-similar nor 
Gaussian (see, e.g., [3). Such models become thus so complicated that they can 
hardly be used to improve the understanding of the underlying phenomenology 
and, at the same time, they are not easily amenable for an analytical treatment. 

The goal of the present work is to give some new phenomenological understanding 
of the problem of turbulent relative dispersion. For this, we make use of large-scale 
direct numerical simulations and essentially focus our study on the role of violent 
events leading to strong fluctuations. To give grounds to possible improvements of 
models for relative dispersion, we introduce some new observables, measure their 
behavior, and interpret the results mostly in terms of geometrical arguments. The 
paper is organized as follows. We first report in Sec. [2] numerical results on the 
time evolution of averaged separations and of higher-order statistics; this is an 
extension of some of our recent work [i^]. In Sec. [3] we focus on the behavior of 
velocity differences and discuss the question of intermittency and of a possible 
multiscaling for the time evolution of their moments. We also introduce a mixed 
moment of separation and velocity difference that attains quickly a scaling regime 
(independent on time and initial separations). In SeclH we apply previous results 
and make use of geometrical considerations to investigate the very large and small 
fluctuations of the inter-trajectory distance. Finally, Sec.[5]encompasses concluding 
remarks and discuss open questions. 



2. Convergence of separation statistics to a scaling regime 
2.1. Numerical results on the mean-square displacement 

We report in this section numerical results on the time evolution of the averaged 
squared distance between tracers. The aim is to provide a first insight on the nu- 
merical simulations that we have used and on the typical values of the timescales 
and lengthscales that are used throughout this paper. One of the major difficulties 
encountered when attacking numerically the problem of turbulent explosive separa- 
tion is that it requires a huge timescale separation. Indeed, to observe Richardson- 
Obukhov superdiffusive regime, one needs to follow particle pairs on a time t 
much longer than the eddy turnover time ~ e^^^^r'^^^ associated to their initial 
separation |-R(0)| = tq and much smaller than the time tl associated to the inte- 
gral scale L. Also, it is in principle required that the initial distance tq belongs to 
the inertial range, so that its associated turnover time has to be much larger 
than the Kolmogorov dissipative timescale = {yje)^!'^. We should thus have 
<^ <^ t <^ tl. To give a concrete idea, if we require that these timescales are 
separated by at least a decade, this implies that one should have a Taylor-based 
Reynolds number R\ = -v/TS (tl/t^) > 4000. Such a high level of turbulence is still 
far from the current state of direct numerical simulations and also from present- 
day experimental setups where accurate particle tracking techniques can be used 
(see [9|, |2]| for recent reviews). For this reason, there is an important need for 



understanding the full process of convergence to the law in order to predict and 
detect for instance possible subleading corrections. This is the spirit in which we 
have tried to present this work. 
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Table 1. Parameters of the numerical simulations 



N 


Rx 


5x 


5t 




V 




e 


^rms 








L 




20483 


460 


3.7-10-3 


6-10-* 


2.5 


10-5 


3.6 


10-3 


0.19 


1.4 


10-3 


0.083 


1.85 


9.9 


4096^ 


730 


1.53-10-3 


1.2-10-3 


1.0 


10-5 


3.8 


10-3 


0.19 


7.2 


10-4 


0.05 


1.85 


9.6 



TV is the number of grid points, Rx the Taylor-based Reynolds number, &x the grid spacing, &t the time 
step, V the kinematic viscosity, e the averaged energy dissipation rate, Mrms the root-mean square velocity, 
r7= (v3/e)-'^/4 the Kolmogorov dissipative scale, r,, = (v/e)-'/^ the associated turnover time, L = «3^-jj,/e the 
integral scale and Tj^ = L/Mrma the associated large-scale turnover time. 



To investigate such issues, we have performed direct numerical simulations of 
the incompressible Navier-Stokes equation in the three-dimensional 27r-periodic 
domain using a standard pseudo-Fourier-spectral solver with a third-order Runge- 
Kutta time marching. Such a method is well adapted to incompressible homoge- 
neous and isotropic turbulence at high Reynolds numbers with an extended inertial 
range of scales. It has the advantages of combining a high degree of accuracy with 
very good performances on massive parallel supercomputers such as BlueGene sys- 
tems and large Intel/ AMD clusters. We have used two sets of simulations whose 
parameters are summarized in table [1] (more details on these simulations can be 
found in [2^). To maintain a statistical steady state, the flow is forced by keeping 
constant the energy content of the two first shells of wavenumbers in Fourier space. 
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Figure 1. (a) Second-order Eulerian structure function S2(x) = (\u{x^r, t) — u{x, t)\^) for the two values 
of the Reynolds number investigated in this paper. The solid line represents Kolmogorov 1941 scaling 
S2(r) ~ (11/3) C2(e»')^/3 with C2 = 2.13. The dashed line corresponds to She-Leveque anomalous scaling 
with ^2 ~ 0.696. The vertical and horizontal colored dashed lines indicate the integral scale and the large- 
scale velocity, respectively, (b) Logarithmic derivative d log 52/d log r. The two horizontal lines encode the 
two scalings shown in Fig. (a). 



The Eulerian second-order structure function measured from these two simulations 
are represented in Fig.[TJ As seen there, the largest one develops a significant scaling 
range where deviations from Kolmogorov 1941 scaling start to be visible. 

For each value of the Reynolds number, the fiow is seeded with 10^ tracer par- 
ticles whose motion is integrated using tri-linear interpolation. After a time suffi- 
ciently long to have converged to a statistical steady state, we start the analysis 
of the dispersion of tracer pairs. For this, we label at a fixed initial time (that 
we fix here to be f = 0) all couples whose distance |-f2(0)| = |Xi(0) — X2(0)| 
is equal to tq ± r/ for tq < 16 ry and equal to vq (1 it 2%) for tq > Wrj with 
ro = 2, 3, 4, 6, 8, 12, 16, 24, 32, 48, 64, 96, 128, and 192r?. This bining was chosen such 
that each family indexed by rg contains a few hundreds of thousands of pairs. We 
then track forward in time all indexed pairs and perform statistics conditioned on 
their initial separation vq. For the sake of simplicity we denote by (•) the Lagrangian 
ensemble average conditioned on tq. 
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Figure 2. Time-evolution of the mean-squared distance for Rx = 730 (a) and Rx = 460 (b) for various 
initial separations ro as labeled. The horizontal and vertical solid lines represent the integral scale L and 
its associated turnover time Tj^, respectively. The dashed line corresponds to the explosive Richardson- 
Obukhov law (jsjl with g = 0.52. 
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Figure 3. (a) Compensated mean-squared displacement {\R{t) — i?(0)p}/(e t^) as a function oit/to with 
to = S2{ro)/{2e) for various initial separations and Rx = 730 (o) and Rx = 460 The two curves show 
behaviors of the form {\R{t) - R{0)\^) ~ get^ + At'^, with A = 52 (ro), given by Batchelor's ballistic 
regime (black dotted line), and A = 2.5/tg (grey dashed line), (b) Measured value of the constant gC in 
front of the sublcading term as a function of the initial separation. It stabilizes to C ~ 1.3/g ~ 2.5 for 
ro ^ »?; this value is represented as a dashed line. 

Figure [2] shows for the two simulations the time evolution of the mean squared 
distance for various values of the initial separation rg. Times and space 

are there represented in dissipative-scale units. After a transient (which roughly 
corresponds to Batchelor's ballistic regime), the mean-squared distance approaches 
the explosive Richardson-Obukhov regime 

{\R{t)f)^get'. (3) 

We observe for both values of the Reynolds number a Richardson-Obukhov con- 
stant g ~ 0.52 ± 0.05. The low accuracy with which this constant is determined 
comes from the fact that, even at the higher resolution, the scaling is observed 
in a rather limited time range. This is even clearer from Fig. [3] (a), which shows 
the compensated mean squared increase of the distance {\R{t) — -R(0)p)/(e i^). On 
this figure, the time has been rescaled by to = 52(ro)/(2e), where ^2 designates 
the second-order Eulerian structure function with absolute values. The choice of 
such a timescale was motivated in [13] as that of deviations from Batchelor's initial 
ballistic regime. Surprisingly, the data shown in Fig. [3] (a) corresponding to var- 
ious initial separations ro far enough in the inertial range seem to collapse. This 
suggests that the timescale to contains most of the dependence of pair dispersion 
upon the initial separation tq. Also data indicate that the subdominant terms in 
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(I3|) are cx t^, leading to postulate g et^{l + C to/t), with a constant C 

independent of ro when ro ^ r]. The product of the constants g and C has been 
estimated numerically and results are shown in Fig. [3] (b) . One can clearly see that 
C w 1.3/5 ~ 2.5 when ro > 8r/. 

One can also see from the figure that C < for rg < 4?]. Thus, for dissipative- 
range initial separations, the asymptotic behavior is attained from below. This 
can lead for such values of ro to an intermediate time range where the mean squared 
distance grows even faster than the explosive law, as for instance observed in 
0]. Another remark that can be drawn from the data is that, independently of the 
Reynolds number, the constant C is equal to zero for ro ~ 4ry. The first subleading 
terms are then oc t, so that the convergence to the law is much faster for such 
an initial separation than for others. This observation could be useful for experi- 
mentalists to optimize their setup. However, such small values of ro are clearly not 
representative of the inertial-range behavior. 



2.2. Higher-order statistics 




We now turn to investigating the large-time behavior of higher-order moments 
of the separation. Figure U] shows the time evolution of {\R{t) — i?(0)|^) (a) and of 
(l-R(t) — R{0)\^) (b). At times smaller than to the separation grows ballistically, so 
that {\R{t) -R{0)\P) ~ tP{\V{0)\P) where V{t) = u{Xi,t) - u{Xi,t) denotes 
the velocity difference between the two tracers. The fact that we have chosen to 
rescale time by to (which depends on second-order statistics of the initial velocity 
difference) implies that the moments do not collapse in this regime because of 
Eulerian multiscaling. However the collapse occurs for t ^ to where these two 
moments grow like t^ and t^, respectively, with possible minute deviations. The 
measured power-laws give evidence that, at sufficiently long times, inter-tracer 
distances follow a scale-invariant law. Also the observed collapses indicate that to 
could be again the time of convergence to such a behavior. 

The presence of a scale-invariant regime is also clear when making use of 
ideas borrowed from extended self-similarity and representing these two moments 
as a function of (|-R(t) — -R(O)p) (see Fig. [5]). This time, for a fixed ro, the 
smallest separations correspond to the ballistic regime. There, we trivially have 
(|i^(t)-l^(0)|P)/(|i^(t)-i^(0)|2)P/2 ~ {\V{0)\P)/{\V{0)\'^)P/^ which has a weak 
dependence on ro, because of an intermittent distribution of Eulerian velocity in- 
crements, but does not depend on time. This normal scaling can be observed for 
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Figure 5. Fourth (a) and sixth (b) order moments of \R{t) — -R{0)| as a function of its second-order 
moment for Rx = 730. The two gray dashed hues show a scale-invariant behavior, i.e. {\R{t) — R{(i)\'^) oc 
(|R(i)-J?(0)p)2 and (\R{t)-R{ld)f) oc (|it(t)-K(0)p)^, respectively. The two insets show the associated 
local slopes, that is the logarithmic derivatives d log(|ii(i) — i?(0)|'')/d log(|ii(i) — i?(0)p), together with 
the normal scalings represented as dashed lines. 

t tQ m. the insets of Fig. O which represent the logarithmic derivatives of the 
high-order moments with respect to the second order. At times of the order of 
noticeable deviations to normal scaling can be observed. Finally, at much larger 
scales, data corresponding to different values of the initial separation rg collapse but 
the curves start to bend down. One observes in the insets that the associated local 
slopes approach values clearly smaller than those corresponding to normal scaling. 
This gives evidence of a rather weak intermittency in the distribution of tracer 
separations. Note that the presented measurements were performed for R\ = 730 
but the same behavior has been observed for Rx = 460. 

To our knowledge, the most convincing observation of an intermittent behavior 
in pair dispersion has been based on an exit-time analysis [23]. However, the rela- 
tion of such fixed-scale statistics to the usual fixed-time measurements we report 
here requires to consider pair separation velocities. As we will see in next Section, 
the velocity difference between two tracers displays statistics that are much more 
intermittent than those for pair separation. This implies that there is no contra- 
diction between an almost normal scaling for distances as a function of time and 
an anomalous behavior of exit times as a function of distance. 
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Figure 6. Probability density function of the distance r at time t = 2.5 to (a) and t = 5to (b) and for 
various values of the initial separation. We have here normalized it by A-kt^ and represented on a log y axis 
as a function of r / {\R(t)\'^)^^'^ . With such a choice, Richardson's diffusive density distribution ^ appears 
as a straight line (represented here as a black dashed line). 



To investigate further this weak intermittency in the separation distribution, we 
have represented in Fig. [6] the probability density function (PDF) of the distance 
for various initial separation and at times where we expect to have almost 
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converged to the explosive regime, namely at t = 2.5 to (a) and t = 5to (b). Such 
measurements are compared to Richardson's diffusive law ([2]). Data suggest that 
a large part of the PDF core (for 0.4 < \R{t)\/{\R{t)\'^)^/^ < 4 at time t = 5to) is 
very well described by Richardson's approach. However, deviations are observed in 
the far tails, at both small and large values of the separation. Such an observation 
is consistent with previous numerical observations P, 23 1. Apparently, these 



deviations affect only weakly the moments we have considered above. We will 
come back to investigating and characterizing them in Sec. |H 



3. Statistics of velocity differences 
3.1. A diffusive behavior? 

In this section, we are interested in the behavior of the velocity difference V{t) = 
u{Xi{t),t) — u{X2{t),t) between two tracers as a function of time. Initially, the 
statistics of V{0) are exactly given by the Eulerian statistics of velocity increments 
at a separation tq. At large times, an explosive behavior for distances implies 
that oc et. A naive picture would consist in interpolating between these 

two behaviors by assuming that 

{\Vit)\^)^S2{ro) + het, (4) 

where h is a positive constant that cannot be straightforwardly expressed as a 
function of As argued in [20] one of the two terms, which correspond to the 
ballistic or to the Richardson-Obukhov regime, is dominant when t is much smaller 
or much longer than to = 5'2(ro)/(2e). As seen from Fig. [7| (a) such a form with 
h = 2.4 (represented as a dashed line) seems to be a good first-order approximation. 
However, this picture is not completely satisfactory. When zooming closer to the 
initial times, one observes that the averaged pair kinetic energy (|^(t)|^) first 
starts by decreasing for a time of the order of a few tens of to- We indeed know 
that initially 

F(t) = V(0)+tA(0) + O(t2), (5) 
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Figure 7. Time behavior of the mean-squared velocity difference for Rx = 730. (a) (|V{t)p), normalized 
by its initial value, as a function of t/to; the gray dashed line is a behavior of the form ~ 
>S'2(''o) + het with h = 2.4. Inset: same but zoomed at short times; the black dotted line is the initial 
decrease associated to kinetic energy dissipation: — ^(fo) — 2et. (b) Time behavior of the 

velocity change - F(O)P); the dashed curve is {\V{t) - V(O)p) ~ S2{ro) + het\ the dotted line has 

slope 2 and corresponds to the initial abrupt variation of the velocity difference. Inset: same but rescaling 
time by r,,. 
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where A(0) denotes the initial difference of acceleration between the two tracers. 
This approximation leads to — »S'2(ro) — 2et (see, e.g., [13] )• According 

to our numerical data, this regime seems to reproduce well what is happening for 
times up to ~ 0.01 to — see inset of Fig. [7| (a). The fact that this prediction ceases 
to describe the data after such a short timescale (that is significantly smaller than 
Iq instead of being of its order) indicates that the subleading terms in ([5]) then 
become important. As a consequence, depending on the quantity we are interested 
in, Batchelor's ballistic regime might end very quickly. 

After this subtle initial decrease, the averaged squared velocity difference ampli- 
tude tends asymptotically to a behavior oc t. However this process occurs in a no- 
ticeably different manner to that of convergence of squared separation to Richard- 
son law: irrespective of the initial separation tq, the quantity (|l^(t)P) always 
approaches from below the asymptotic behavior. In other terms, we observe that 
for a fixed value of t/to, the pair kinetic energy (| V(t)p) is an increasing function 
of ro/r? and seems to converge to (jl]) only in the limit ro/rj — )• oo. 

At large times, the mean-squared velocity difference grows linearly. This be- 
havior, together with the fact that acceleration differences are correlated over 
timescales of the order of the Kolmogorov dissipative time r^, suggest that ve- 
locity differences have a diffusive behavior for t ^ r^j (see [i^]). If this was exact, 
it would lead to the behavior (j4]) as the temporal increments of V would be inde- 
pendent of its initial value. Also, a purely diffusive behavior of V would result in 
the fact that its mean-squared temporal increment {\V{t) — V{0)\'^) is (x t for all 
times t ^ Tfj. As seen in Fig. [7](b), the diffusive regime is actually only asymptotic. 
Also, data show that V{t) initially changes on times of the order of by a factor 
of the order of its typical initial value [5*2 (ro)]^/^. This abrupt evolution can be in- 
terpreted phenomenologically. With some finite probability, one of the two tracers 
is within a vortex filament at time to- The typical energy content of this filament 
will contribute to the value of 52(ro). However, after a time t of the order of r^, the 
trajectory of this tracer will have turned around this filament, so that its velocity 
will have completely changed orientation (without changing much its amplitude). 
This will result in |V(t) — ^(0)| ~ |^(t)|, explaining the observed behavior. As a 
result of this sudden kinematic variations of velocity differences, {\V{t) — V{0)\'^) 
behaves in a very similar manner to (|^(t)P) for times t ^ (compare the dashed 
lines in Fig. [7|(a) and (b)). Let us also notice that the convergence to this behavior 
is again from below, irrespective of the initial separation rg. 

There is hence an abrupt change (occurring on timescales of the order of t.^) 
that prevents from determining an effective initial velocity difference and thus 
from observing a clear diffusive behavior of V{t). However, data suggest that the 
timescale of convergence to this behavior is, as for separations, of the order of to- 
To understand further this question, we next turn to investigating the behavior of 
the longitudinal velocity difference between the tracers. 

3.2. Geometry of longitudinal velocities 

We are here interested in the evolution of the longitudinal component V^"(t) = 
R{t) ■ V{t)/\R{t)\ of the velocity difference as a function of time. This quantity is 
important to characterize pair separation as d|i?|/dt = V^^. Initially, the averaged 
longitudinal velocity vanishes, i.e. (1^11(0)) = 0; this is due to the statistical sta- 
tionarity of the fluid flow. For times t ^ to in the Batchelor's ballistic regime, the 
pairs are keeping their initial velocity difference and one can easily check that 



{vHt))c^t{\V^{t)\')/ro, 



(6) 
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Figure 8. (a) Time evolution of the averaged longitudinal velocity difference (V'l(t)) for various initial 
separations, as labeled, and for = 730. (b) Averaged time t» = —rovo/(\V{0)\'^ 11/11(0) = vq) when 
trajectories with a given initial longitudinal velocity difference vg reach a minimal distance. 



where V denotes the components of V that are transverse to R. It is thus clear 
that the average velocity at which tracer trajectories separate immediately be- 
comes positive. Figure [8] (a) , which represents the time evolution of (V^ll(i)), shows 
without doubt this initial linear growth. This increase of the longitudinal velocity 
difference has an interesting geometrical interpretation. If all pairs were to main- 
tain indefinitely their initial velocity differences V{0), it is clear that they would 
eventually all separate at large times. Indeed, they will reach a minimal distance 
at a finite time = -[R{0)-V{0)]/\V{0)\'^ = -rQV\\{0)/\V{0)\'^ , which is positive 
for particles that are initially approaching. After the time t*, the distance between 
particles increases and V^^ becomes positive. This leads to an increase of (yll) that 
comes from kinematic considerations and is obviously not due to any dynamics im- 
posed by the turbulent flow. To estimate the typical value of the minimal distance 
at time t*, we have performed statistics on pair separation conditioned not only 
on the initial distance rg, but also on the initial longitudinal velocity difference 
by binning pairs with 1/11(0) = ± 6vq. From these statistics we have defined the 
averaged time t^: = — rofo/(| V(0)p \vo) at which trajectories with a given initial 
longitudinal velocity difference vq are at a minimal distance. The data are shown 
in Fig. [8] (b) for various values of tq and as a function of vq. We observe that for 
ro ^ rj this time approximatively takes the form t^, ~ to f{vo/ {\V {0)\'^)^^'^), where 
the function f{x) attains its maximum (roughly equal to 0.05) at x ~ 1. This 
confirms the observation made in Fig. [8] (a) that the initial growth of (V^") occurs 
on a time length of the order of a few hundredths of to. 

In addition to the change in the mean longitudinal velocity difference discussed 
above, numerical measurements show that the full distribution of looses its 
symmetry and develops fatter tails when time increases. Figures [9] (a) and (b) 
represent the skewness S and flatness of as a function of time and for the 
same initial separations as in Fig. [8] (a). These observables are frequently used in 
turbulence to quantify the shape of the velocity increment distribution. For two- 
particle Lagrangian statistics, they are defined as 

(ivHt)-iv'mi (y"w-(v"(W 

<iv"(«)-(v"(tm-" ([v"(t)-(v'(<))P) 

As already observed for instance in [24] , these quantities strongly vary as a function 
of time and maintain a marked dependence upon the initial separation ro for rather 
long times. Figure [9] (a) shows that the skewness of V^"(t) starts from negative 
values (to be in agreement with the Eulerian 4/5 law) and becomes positive at 
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Figure 9. Time evolution of (a) the skewness S{t) and (b) the flatness J^{t) deflned in JTJ for the same 
initial separations ro as in Fig.|8](a). The dashed lines are the values corresponding to a Gaussian distribu- 
tion: 5 = and = 3. The dotted line in (a) is the initial value of the skewness obtained when assuming 

3/2 

Kolmogorov 1941 scaling, namely S = {4/5)/C2 . 

times larger than « 0.01 to- This initial change of sign can also be interpreted 
geometrically in terms of the time when initially ballistically approaching pairs 
begin to move away. However, after this, the curves separate and each of them 
attains a maximum at times of the order of to or slightly smaller. This maximal 
value of the skewness strongly depends on the initial separation: less is ro, higher 
it is. After this maximum, the skewness decreases without attaining an asymptotic 
regime that would be independent of tq. This could be either due to the fact 
that there is a persistent memory of ro in such quantities or to a contamination 
by finite Reynolds number (and finite size) effects. The same kind of behavior is 
observed for the flatness of the distribution of as seen in Fig. [9] (b). However, 
the increase of its maximal value when decreasing ro is even more pronounced. 
Another noticeable difference is that the initial value of J- itself depends on ro and 
relates to the scale dependence of the Eulerian flatness. 




10""' 10 " 10 ' 10° lo' 10 10 ^ \ 10"' 10° lo' 



Figure 10. Time evolution for ro = 24r; of the averaged distance (a) and of the squared velocity difference 
(b) conditioned on yll(O) = do with bo (0)]2)l/2 = 0.1,0.4,0.8, 1.2, ...,2.8. The red curves correspond 
to positive values of vq and the blue ones to negative values. In (b) the dashed-dotted line is ex f 

The strong dependence on ro of the skewness and of the flatness can be inter- 
preted phenomenologically in terms of the intermittent nature of velocity incre- 
ments. A large positive value of S corresponds to a large probability of having 
pairs separating faster than the average. Such events will also be responsible for a 
large value of the flatness It is clear that those particles that separate the faster 
are typically those which are the most separated. Also, in a turbulent flow, the 
typical value of the velocity increases with scale, so that the particles which get 
quickly separated are likely to continue separating faster than the average. This is 
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evidenced in Fig. [TOl which represents the mean squared distance and the averaged 
squared amphtude of the velocity difference conditioned on the initial value of the 
longitudinal velocity difference for ro = 24?;. One observes that there is up to the 
latest time (of the order of 5 to); a noticeable memory on the initial value of 
Besides this consideration, one also remarks in Fig. [10] (b) that the pairs having an 
initially large negative longitudinal velocity difference (the blue curves) dissipate 
much more kinetic energy than the others. However, this does not prevent them 
from separating at large times faster than the pairs having initially a smaller veloc- 
ity difference. We now turn back to the explanation of the long-term dependence 
on ro of the skewness and flatness of Recall that, in a turbulent flow, violent 
velocity differences are more probable at small scales than at larger scales. This 
implies that pairs with a small initial separation are more likely to experience a 
large (positive or negative) initial velocity difference. This will make them separate 
faster than the average and thus experience even larger values of the velocity. The 
rapid and strong increase of fluctuations in their velocity differences is thus due 
to a kind of snowball effect. We will come back later to quantifying with more 
accuracy this phenomenon. 

3.3. Intermittency 

All the considerations on the dependence of V^^ upon the initial separation are also 
visible in the probability distribution of the longitudinal velocity. Figure [TT] (a) 
shows the centered PDF of V^^ normalized to unit variance for various times and 
ro = 12ry. The data clearly show that at times later than 0.01 to, there is a change 
in the sign of the skewness. Also, one sees that the time dependence of the skewness 
and of the flatness comes from the right tails associated to large velocity differences, 
supporting the arguments discussed above. The left tails, which corresponds to 
approaching pairs, seem on the contrary to collapse. The right tail has a very rich 
behavior. It starts with broadening at times t < to, in agreement with the increase 
of flatness. For t > to, it then decreases and possibly goes back asymptotically to 
its initial form. If this was true, it would imply that the distribution of velocity 
differences keep in memory the initial separation at any later time. 




[V"(t) - <v"(t)>] / <[v"(t) - <v"(t)>]->"^ [V"(t) - <v"(t)>] / <[v"(t) - <v"(t)>]^>' 



Figure 11. Centered and normalized probability density functions of the longitudinal velocity difference: 
(a) for ro = 12?7 and various times and (b) for various initial separations and t = 5to. In both cases, the 
black dashed line shows a Gaussian distribution. 

This is also evidenced from Fig. [11] (b) , which shows the same PDFs for various 
separations and a time t = 5tQ fixed. Again we observe a rather good collapse of 
the tails associated to negative longitudinal velocity differences, but the right tails 
display very strong dependence on the initial separation. Clearly, the behavior of 
this tail is a stretched exponential for ro ^ Srj and is faster than exponential for 
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< r|j = 64Tl 



Figure 12. Fourth-order (a) and sixth-order (b) moments of the longitudinal velocity difference as a func- 
tion of its second-order moment for various times and initial separations. The two dashed lines correspond 
to a scaling compatible with that of Lagrangian structure functions proposed in 1251 , namely C4 /C2 ~ ^-71 
and C6/C2 = 2-16. The insets show the logarithmic derivative d log([Vll (t)]J')/d log{[y H (t)]^) for (a) p = 4 
and (b) p = 6 as a function of t/to; there the bold dashed lines show the Lagrangian multifractal scaling 
and the thin lines what is expected from a self-similar behavior. 



larger initial separations. The actual level of statistics do not allow us to relate 
systematically this behavior to that of the initial velocity difference distribution. 

Finally another way to address the question of intermittency of the velocity 
difference consists in finding how moments of its longitudinal component depend 
on time. For that we follow, as in the case of the moments of distances, an approach 
similar to that of extended self-similarity. Figures [T2] (a) and (b) show the fourth 
and sixth-order moments of as a function of its second-order moment. As 

evidenced in the insets, they display an anomalous behavior that differs from simple 
scaling. However, the collapse for various tq is much less evident than for the 
moments of the distance, except perhaps at sufficiently large times. One can there 
guess a power-law dependence of ([^"(i)]'') as a function of Surprisingly 
the power is compatible with the scaling exponent of the Lagrangian structure 
functions that were obtained in [isl] by relating velocity increment along trajectories 
to She-Leveque multifractal spectrum for the Eulerian field. The two dashed lines 
in Fig. [12] (a) and (b) correspond to the predicted values C4 /C2 — l-'''^ ^^'^ Ce /C2 — 
2.16. Confirming further this match would require much better statistics. 



3.4. Stationary distribution of rescaled velocity differences 

As we have seen previously, the velocity difference between tracers displays very 
intermittent features and, as a consequence, does not converge to a behavior with 
temporal self-similarity, or does it only very slowly. The situation is very different 
when interested in mixed statistics between distances and longitudinal velocity 
differences. As seen in [2^, the moment ([y" (t)]^/|-R(t)|), which is initially negative 
and equal to — (4/5)e, tends very quickly to a positive constant — see Fig. [13] (a). 
The decrease at very large times comes from the contamination of the statistics by 
pairs that have reached a distance of the order of the integral scale. The asymptotic 
value ~ 6.2 e seems to depend only weakly on the Reynolds number. The collapse 
of the curves associated to different Reynolds numbers and, for tq ^ ry, to various 
initial separations indicate that the time of convergence is oc to- Figure [13] (b) 
shows the same moment but conditioned on the sign of the initial longitudinal 
velocity difference. One observes that for initially separating pairs (red curve), the 
convergence to the asymptotic value is on a time of the order of r^. Conversely 
for tracers that initially approach each other (blue curve), the convergence is less 
fast. We have seen in Sec. 13.21 that such pairs first attain their minimal distance at 
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Figure 13. (a) Time evolution of the mixed moment ([V" (t)]^ /\R{t)\) for different initial separations and 
the two values of the Reynolds number: Rx = 730 (o) and Rx = 460 {+). The black dashed line shows 
the asymptotic value 6.2 e. (b) Same for ro = 24r; and Rx = 730 in black and for the same moment but 
conditioned on positive (red curve) and negative (blue curve) values of the initial longitudinal velocity 
difference ^11(0). 

tKit^^ 0.01 to- Then, at that time, {[V^t)f /\R{t)\ \ V^O) < 0) changes sign and 
the convergence to ~ 6.2e occurs only later. Such initially approaching pairs are 
leading the average so that the overall convergence is on timescales of the order 
of t^,. The mixed moment which is a kind of "local dissipation" along 

pairs of trajectories, is thus conserved by the Lagrangian flow. 

Actually, it is not only the average of the "local dissipation" that converges to a 
constant but its full distribution seems to attain a stationary regime on times of the 
order of t^. Figure [14] (a) shows for a given initial separation, the convergence at 
large times of the PDF of [yll(t)]'^/|/?(t)|. One observes that the right and left tails 
converge on different timescales. The tail associated to large positive values (sepa- 
rating pairs) occurs on timescales of the order of (which is in this case ~ 0.07 to), 
while that for negative values (approaching pairs) converges slower. For the largest 
time {t ~ 5 to) one observes that deviations to the asymptotic distribution occur 
again at very large positive values. This is due to a contamination of such events 
by the large scales of the turbulent flow. This decrease is in agreement with the ob- 
served departure in Fig. [13] of the average from its asymptotic value at large times. 
Figure [HI (b) shows the PDFs of for different initial separations and at 

a fixed time sufficiently large to be ensured that all distributions have attained 
their asymptotic regime. One observes a robust collapse, much more pronounced 
than for both the distribution of separations and that of velocity differences. Note 
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Figure 14. (a) PDF of the "local dissipation" [yH (t)]"^/|ii(t)| for ro = 24r; and various times (as labeled), 
(b) Same for various initial separations and a time t = 20 Tri fixed, t/tg goes here from 0.2 to 5 so that, 
in all cases, the mixed quantity has reached its asymptotic regime. In both figures the vertical 

dashed lines show the position of the average value fs 6.2e. 
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that in Fig. [T3] the distributions are raw and were not rescaled by any moment of 
The asymptotic PDF is peaked around zero (rather than its mean value) 
and displays fat tails that, according to our data, are oc exp(— C on 
both sides. 

To our knowledge, such a fast and manifest convergence of these mixed statistics 
has never been reported before. From a phenomenological viewpoint, one expects 
oc t oc in the explosive Richardson-Obukhov regime, so that the "local 

dissipation" should become constant at sufficiently large times. However 

we have observed here that the convergence of this quantity to its asymptotic value 
occurs much faster and in a much more definite manner than the convergence of the 
statistics of \R\ and |V| to their respective asymptotic forms. This indicates that 
the statistical stationarity of the "local dissipation" is more likely to be a cause 
rather than a consequence of Richardson-Obukhov explosive separation. At the 
moment we unfortunately lack a clear understanding of the physical mechanisms 
that are responsible for the observed behavior of and that could shed 

light on its relation to explosive separation. 



4. Extreme events in separation statistics 

In this section we turn back to the statistics of the distance between tracers. 

Our goal is to explain the mechanisms leading to very large or very small values 
of this distance in the light of the various observations made in previous sections. 
Considering the relative dispersion of tracers conditioned on their initial distance 
tq can be geometrically interpreted in terms of the time evolution of an initially 
spherical surface of radius ro that is centered on a reference trajectory. The trans- 
port of this surface by a turbulent flow is generally very complex. Incompressibility 
implies that the volume of the sphere is conserved but the velocity field roughness 
will be responsible for strong distorsions of its surface. This is represented in Fig. [15] 
that shows for Rx = 460 at three various times the shape defined by the instan- 
taneous position of 60 trajectories that are all initially at a distance 24r/ from a 
reference tracer. We can qualitatively deduce from these snapshots that the large 




Figure 15. Quasi-Lagrangian evolution of the set of tracers that are initially at a distance ro = 24r; of a 
reference trajectory. 

excursions of the inter-trajectory distance go together with strong pinches of the 
surface. We also see that the intense stretchings are occurring in a time-correlated 
manner: the surface angles that are visible at late time have formed at very early 
stages. As we will now see, the most-separated pairs have been so for long times 
and carry a reminiscent dependence on the initial separation. 
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4.1. Memory in large- distance statistics 

We first consider the events related to distances that are much larger than their 
average. It was argued in Sec. 12.21 that the large-value tail of the separation PDF 
is not well described by Richardson's distribution. As seen in Fig. [6] for t = 5to; 
the tail is indeed broader than exp(— C for tq < 8r] and narrower otherwise. 

It seems to tend to a Gaussian when either rg or t are sufficiently large. We have 
checked that these behaviors are not due to a contamination by pairs that have 
reached the large scales of the turbulent flow, as the considered distances are still 
well below L, except maybe for the largest initial separation tq = 24t/. These 
observations suggest that for such extreme events, there is a long-term memory of 
the initial separation. 

To qualify further the history of pairs that are well separated at large times, 
we have performed the following analysis. Fixing a final time tf sufficiently large 
to have reached the explosive Richardson-Obukhov regime, we have carried out 
statistics conditioned on pairs that are far separated at t = tf, say such that 
their distance is |i?(tf)| > 2{\R{t{)\). In order to not be contaminated by finite 
inertial subrange effects, we have restricted this analysis to Rx = 730 and to initial 
separations 2 7] < < 24ry, and we have chosen the largest compatible value of tf, 
namely tf = 5 to- Let us denote by (•)+ the resulting conditional ensemble average, 
i.e. (•)+ = (• I {|i?(tf)| > {\R{t{)\)}) . Figure [TBI (a) represents the relative increase 
(|i?(t) — ii(0)|^)+/(|i?(t) — i?(0)p) of the mean-squared displacement. The various 
curves, which are associated to different initial separations tq, have a maximum at 
t = tf = 5 to represented by a dashed line. The value of this maximum is becoming 
larger when tq decreases and approaches the dissipative scales. This is a signature 
of the dependence of the large-positive value tail upon r^. On the right-hand side 
of the maximum, the various curves tend back to 1. This indicates that separations 
that are large at a given intermediate time relax in average to an ordinary behavior 
at later times. This can be regarded as a consequence of the weakening of large- 
distance probability tails as a function of time. The situation is rather different 
when interested in the left-hand side of the maximum. The curves do not converge 
to 1, and thus to an average behavior when t — )• 0. This means that pairs that are 
well separated at a given time are likely to have been so at any previous times. 
This asymmetry when going forward or backward in time could already be grasped 
in Fig. [T5j The trajectories that are far away from the reference tracer at the latest 




Figure 16. Conditional statistics over pairs that are at a distance at least twice the average at time 
tf = 5to (conditional average on |fi(5to)| > 2(|ii(5to)|) is denoted by (■) + ). (a) Relative increase in 
the mean squared change of separation (\R{t) — J?(0)p) as a function of time, (b) Averaged longitudinal 
velocity difference conditioned on the separation at time tf = 5to (+) and without conditioning (solid line). 
Inset: initial value of the conditioned longitudinal velocity difference as a function of the initial separation; 
the solide line is oc r„ . 
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time (as, e.g., those defining the left and top corners) are visibly also well separated 
at the intermediate time. Conversely the right corner at the intermediate time have 
later stopped separating faster than the average. 

Figure [16] (b) shows the time evolution of the averaged longitudinal velocity 
difference conditioned on having a separation twice the average at time 

tf = 5 to- A first remark that cannot be seen from the log- log plot is that the 
initial value of this average is strictly positive. The inset represents the variation 
of (yl'(0))-|- as a function of the initial separation. This implies that the pairs 
that are well separated at a late time are preferentially separating from the very 
beginning. The initial longitudinal velocity fluctuation that is necessary for the 
pairs to be far apart at time tf becomes weaker when ro increases. Data suggest 
that (yll(0))+/[5'2(ro)]"^/^ oc Tq '^^^ as seen in the inset. This initial separation makes 
the selected pairs reach an almost diffusive regime at time much shorter than the 
end of the average ballistic regime. For large-enough initial separations we indeed 
have (yll(t))+ ~ for O.Uq < * < ^o- The pair distance encountered a final 
acceleration right before the time of conditioning tf. After that, the longitudinal 
velocity difference relaxes slowly to the average regime. Again here, as in the case 
of the average separation, we observe that the imposing of having a large distance 
at a large time selects pair histories. The main contribution to statistics is indeed 
given by couples encountering an initially violent separation and quickly reaching 
larger scales. Such pairs converge rapidly to an explosive regime and continue to 
separate faster than the average for a long time. This is the snowball effect we have 
mentioned before. 
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Figure 17. Exponent a of the right tail of the PDF of distances as a function of time and for 

various initial separations. The exponent was obtained by fitting — logp(|i?l) to a power-law for {\R{t)\) < 
\R{t)\ < 0.5 L. 

In order to give a more quantitative handle on the far tail of the separation PDF, 
we have estimated and fitted its functional shape and studied its variations. For 
that, following the observations made in Sec. 12.21 we have assumed that p(|-R|) oc 
exp(— C |-R|") and measured how the exponent a depends on both time and initial 
separation. Figure [T71 shows the time behavior of the exponent a for various values 
of the initial separation. On the figure, the three dashed horizontal lines represent 
the tails of a Gaussian (a = 2), of an exponential (a = 1) and of Richardson's 
distribution [a = 2/3). At time t = 0, the distribution is peaked around ro and 
compactly supported, so that a = oo. The exponent then decreases, crosses a = 1 
(so that the distribution becomes a stretched exponential) for t ~ 0.1 to, and reaches 
a minimal value that depends on the initial separation ro. For the smallest ro, 
this minimal value is below Richardson's prediction, as previously noticed. For 
larger separations, this minimum increases and when vq ^ rj the curves seem to 
saturate to the value a = 1. Also, we cannot exclude that all curves converge to 
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the exponential value when t — )• oo. The increase at the last stage can hardly be 
blamed on an integral-scale effect. We have indeed excluded here all pairs that are 
separated by a distance larger than L/2. 



4.2. "Fractal distribution" at small distances 




Much more violent and intermittent events take place for particle pairs that sep- 
arate much less than the average. As already stressed and observed in Fig. [T5l the 
large excursions of inter-trajectory distances go together with strong pinches of sep- 
arations. This is evidenced from Fig. [18] that represents the probability distribution 
of |-R| at various times and for two inertial-range values of the initial separation. 
One observes that at very short times, the PDF is peaked around tq. When time in- 
creases, its maximum, which roughly corresponds to the value of shifts to 
larger values and the distribution broadens simultaneously at large and small val- 
ues. This leads to the development for ro <C |-R| <C (|-R|) of an intermediate range 
of pairs whose separations lag behind the average evolution. In this subrange, the 
PDF behaves as a power law oc |-R|^, where the exponent /3 evolves as a 
function of time and tq. The power-law behavior is substantiated when measuring 
the cumulative probability P^{r) = Prob[|i?| < r] = p[r')dr' oc r^~^^ of inter- 
trajectory distances. One expects from Richardson arguments that /3 = 2, so that 
P'^{r) oc r^. 

The exponent /3 can be interpreted in terms of fractal geometry. If all the tra- 
jectories were uniformly and independently distributed in space, the fraction of 
pairs at a distance less than r would be oc r'^, where d = 3 is the space dimen- 
sion. Also one would expect this fraction to be oc r'^~^ if the trajectories were all 
confined on a surface, and oc r'^"^ if they were on a curve. In general, the expo- 
nent /3 + 1 = limr^.o[log-P^(?')]/[logr] measures the correlation dimension of a 
fractal set. Richardson's behavior thus corresponds to the idea that, at suffi- 
ciently long times, trajectories forget about their initial separation and distribute 




Figure 19. Two dimensional sketch of the time evolution of trajectories that are initially located at a 
distance ro of a reference tracer (shown as a black dot at the center of the circle) and then spread. 
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homogeneously in space. As we have previously discussed, pair dispersion can be 
geometrically interpreted as the average Lagrangian evolution of an initially spher- 
ical surface of radius ro centered on a reference trajectory. At time t = 0, it is 
clear that P^(r) is a Heaviside function centered on r = rg, so that /3 = oo. At 
later times the stretching and pinching of the sphere leads to a non-trivial be- 
havior of P^(r) for r <^ {\R\)- This behavior is sketched for two dimensions in 
Fig. [T9l The exponent /3 is thus measuring the fractality of the image of the sphere 
by the Lagrangian flow. It relates to the fraction of this forward-in-time image 
that remains within a distance r of the reference trajectory. At sufficiently large 
times (right-most panel of Fig. I19p . one expects the set of trajectories to be mainly 
stretched by the flow eddies with the most relevant correlation time, that is those 
of size (l-R(i)l), and the small scale fluctuations to become less and less important. 
Hence, at very large times, the only pairs that are still at a distance r <^ {\R{t)\) 
should be distributed on a surface (curve in two dimensions), so that /3 -|- 1 ~ 2 and 
p(r) oc r. This approach contradicts the prediction of Richardson that, somehow, 
postulates that the velocity difference correlation time does not depend on scale 
(as it is assumed to be zero). 





Figure 20. (a) Cumulative probability distribution of the distance P^(r) = Prob[|i?| < r] as a function 
of r for ro = 4»7 and at various times; the gray dashed curve is the cumulative probability associated 
to Richardson's distribution that displays a behavior oc at small distances, (b) Time evolution of the 
left-tail power-law exponent /3 for various values of the initial separation. The black solid line shows /3 = 2 
corresponding to Richardson's distribution and the dashed line stands for the value /3 = 1 that is argued 
in the text. 



To shed light on such ideas, we have estimated the exponent /3 for various ro 
and at different times. For this purpose, we have used the cumulative probability 
P^(r) (shown in Fig. [20] (a) for rg = 4r/) to obtain P + 1 through a fit of its 
local slope. The resulting measurements are displayed in Fig. [20] (b) as a function 
of time. Up to the large error bars due to a lack of statistics, one observes that 
the timescale to is again relevant to describe both the initial dynamics and the 
convergence to an asymptotic regime. At small times t <^ to, the exponent /3 takes 
large positive values, which, as explained above, are due to the initial conditioning 
|i?(t)| = ro- At time t ~ 0.3to) the exponent becomes smaller than 3. For t ^ to, it 
asymptotically approaches the value /3 = 2, so that data give evidence in favor of 
the argument exposed above. Interestingly, for small initial separations tq < 16r/, 
the limit is reached from below and preceded by a minimum of /3. This means 
that for such intermediate times, the pinching (or equivalently the stretching) of 
the sphere is so strong that the intersection of the latter with spheres of radius r 
defines a fractal object that is more concentrated than a surface. Such a temporary 
singular behavior is certainly related to the presence of tails fatter than exponential 
at very large distances. 
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5. Concluding remarks 

In this paper we have confirmed and extended previous results of [20] on the 
timescales of convergence of turbulent pair dispersion to an asymptotic regime. 
We have seen that low and medium-order moments of the particle separation ap- 
proach an asymptotic regime on times of the order of to = S2{ro)/ (2e), where 52(ro) 
denotes the (absolute value) structure function associated to the initial separation 
ro and e is the mean rate of kinetic energy dissipation. This timescale has been 
shown to be relevant to describe also the initial kinematic change of sign of the lon- 
gitudinal velocity difference and the tails of the distance distribution at small 
values. However, we have seen that to is not relevant to describe the convergence of 
the velocity difference statistics to an asymptotic regime. We have observed that 
up to the largest time, the skewness, the flatness, and more generally the shape of 
the velocity difference distribution still depend on the initial separation, even when 
time is rescaled by to- This leads to a behavior that is by far more intermittent 
than that of separations. Also we obtained evidence that the far tail distribution 
of separation is also keeping a nontrivial memory on ro up to the largest times. 
Throughout this paper, we made an important use of geometrical considerations to 
explain phenomenologically the statistical events leading to extreme fluctuations. 
In particular, we argued and obtained numerical evidence that, at sufficiently large 
times, the probability distribution p(r) of inter-tracer distances is oc r for r ^ 
and decays as exp(— Cr) when r 3> These two observations strongly contra- 

dicts Richardson's eddy diffusivity approach, which seems nevertheless to give a 
good approximation of the core of the distribution. Finally, we obtained a striking 
result concerning mixed distance and velocity statistics. We indeed found that the 
distribution of a "local dissipation", defined as attains relatively quickly 

a scaling regime that is independent of both tq and time. This behavior, which, 
to our knowledge, has never been observed before, gives very strong constraints 
for the development and validation of stochastic Markovian models for turbulent 
relative dispersion. 

A central question that we raised concerns the physical mechanisms leading to 
this fast convergence of the "local dissipation" to a statistically stationary behavior. 
At the moment, an even phenomenological explanation is still missing. We have 
understood that kinematic considerations can be used to explain why 
which is initially negative, becomes quickly positive. This is due to the fact that 
initially approaching trajectories will always eventually separate. However, there 
is no clear idea of why this average should converge to a time-independent value. 
Possible ideas could consist in finding other mixed moments that are strictly (and 
not asymptotically) conserved by the two-point turbulent Lagrangian flow. Such 
considerations are kept for future work. Finally, another possible extension of the 
current study is to apply the developed understanding of the geometry of relative 
dispersion to more complicated turbulent transport situations involving more than 
two trajectories. This is for instance the case for the forward-in-time dynamics 
of triangles or tetrahedrons studied in [2^. The quality of today numerical data 
would be very useful to revisit such questions with an emphasize on extreme events. 
Another situation that is much closer to applications is that of the transport by 
turbulence of pollutant patches. We have seen in the current study that relative 
dispersion relates to the distorsion of a sphere by the Lagrangian flow. We related 
the small-scale behavior of separation statistics to the fractal dimension of this 
object. Such geometrical considerations could be generalized to understand and 
quantify the large concentration fluctuations in the dispersion of a pollutant spot. 
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